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We present an effective medium theory in order to ana- 
lyze the reported optical properties of aligned carbon nan- 
otube fihns (W.A. de Heer et al, Science 268, 845 (1995)). 
This methodology is based on photonic band structure cal- 
culations and allows treatment of complex media consisting 
of particles that interact strongly. We also develop a simple 
Maxwell- Garnett type approach for studying this system. In 
comparing the results of both mean field theories, we demon- 
strate that the inclusion of the full electromagnetic coupling 
between the nanotubes, as our numerical scheme does, is nec- 
essary for a complete explanation of the experimental data. 
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Since the discovery of tubular fuUerenes in 1991 [y, 
several experimental and theoretical studies of the unique 
properties of these materials have been carried out. De- 
spite the short time since the first carbon nanotube was 
synthesised, possible applications of these new materials 
have already been reported. Nanotubes can be used as 
atomic-scale field emitters |^ or as a pinning material in 
high-Tc superconductors |^] whereas it has been proposed 
that hollow tubes might serve as nanoscale molds. The 
unique properties of these carbon tubes are derived from 
their special topology and their nanometric dimensions. 
Very recent measures of the conductivity of individual 
nanotubes suggest that the geometry of these tubes 
plays a fundamental role in determining if the nanotube 
is metallic or semiconducting. From the mechanical point 
of view, nanotubes seem to be very strong and flexible at 
the same time due to the ability of sp^ carbon atoms to 
rehybridize when the bonds are deformed. 

These carbon nanotubes also present very interesting 
optical properties. Thin films of aligned carbon tubes are 
birefringent, reflecting differences in the dielectric func- 
tion for light polarized along (s-polarization) and normal 
(p-polarization) to the tubes g]. Up to our knowledge, a 
complete theoretical analysis of the optical properties of 
these arrays of carbon nanotubes is still missing. 

In analyzing the propagation of electromagnetic waves 
in heterogeneous media, various effective medium ap- 
proaches like the Maxwell-Garnett (MG) approximation 
1^ have been used to determine the effective dielectric 
constant of different composite structures. Very recently, 
the validity of the MG approach has been checked by 
comparing their results for periodic dielectric systems 
with the ones obtained in the long-wavelength limit of 



photonic band structure calculations ||^ . For these struc- 
tures, the MG results seem to agree reasonably well with 
their vector- wave counterparts for small flUing ratios, but 
in the close-packing regime the inclusion of the full elec- 
tromagnetic (EM) interaction between the different ob- 
jects is unavoidable. Only numerical methods give the 
capacity to treat this more complex limit. 

Our idea is to apply a numerical methodology orig- 
inally constructed to study photonic materials and the 
propagation of EM waves in complex media to an- 
alyze the optical properties of these carbon nanotube 
films where the tubes form a very close-packed structure. 
Within this framework it is possible to develop an effec- 
tive medium theory in which the EM interaction between 
the constituents of the complex system is fully included 
Q . We also want to develop a MG-like approach able to 
study arrays of carbon nanotubes and compare these re- 
sults with our exact results, studying the range of validity 
of MG approximation for carbon based materials. 

Our model for a carbon nanotube film is shown in 
figure la. An array of tubes, infinitely long in the z- 
direction, is arranged on a square lattice. The diameter 
of the cylinders is chosen to be 10 nm, in accordance 
with experimental evidence Q. The distance between 
the nanotubes, d, will be varied throughout our calcula- 
tions, studying the dependence of the optical properties 
of this structure on the separation between the tubes. 

Let us discuss first the dielectric model we are using 
for describing the dielectric characteristics of an isolated 
carbon nanotube. Garbon nanotubes are made of a num- 
ber of cylindrical shells of planar graphite. Graphite is 
highly anisotropic and its optical response depends on 
the direction of the electric field with respect to the axis 
normal to the graphite sheets. Therefore, for a given 
frequency it is necessary to distinguish two different 
components, e||(w) and e^{uj), in its dielectric function 
for the directions parallel and normal to the axis, respec- 
tively. These components have been measured by several 
methods, leading to slightly different sets of data. For 
our calculations we have used the dielectric functions of 
graphite as tabulated in Ref. ||l^. In transfering these 
dielectric constants to cylindrical multishells we will as- 
sume that this material is a dielectric continuum and lo- 
cally identical to graphite. A similar transfer procedure 
has been previously applied to compute the UV absorp- 
tion of multishell fuUerenes ||ll[. At every point inside 
the nanotube we can write down a local dielectric ten- 
sor that in cylindrical coordinates (see fig lb) takes the 



diagonal form: 



e((jj) = ej_{iLj){69 + zz) + e|| ((jj)rr. 
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In Cartesian coordinates, the tensor like character of 
the dielectric function is more pronounced: 
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Arc-generated multilayered nanotubes are generally 
hollow. Their inner radius distribution usually varies 
from 0.25 to 2 nm, being tubes with inner radius of 0.75 
to 1 nm the most abundant ones in the samples p2[ . 
The fundamental importance of the hollow character of 
isolated nanotubes in their dielectric response has been 
previously stressed [|l3|. However we have found that 
when nanotubes form close-packed structures, as in nan- 
otube films, the presence of a hollow core plays a much 
less important role. Moreover, our results indicate that 
in the strong coupling limit and for p-polarized light, the 
hollow core can be safely ignored if the ratio between the 
inner and the outer radii of the tubes is less than 0.4. As 
the nanotubes present in the experiment |^ seem to ful- 
fill this condition, for p-polarized light we will only show 
results for plain cylinders (inner radius equal to zero) . A 
detailed discussion of the effect of a hollow core on the 
dielectric properties of aligned carbon nanotubes will be 
presented elsewhere 

If the nanostructured film could be replaced by an 
effective homogeneous medium, a single pair of Bloch 
waves for each polarization, {ks{Lu),—ks{uj)) and (kp{Lu),- 
kp{u})), should dominate the photonic band structure of 
the system. It is indeed the case that over the entire 
range of frequencies w, four EM waves are dominant and 
we can use the dispersion relation of these Bloch waves 
to define the effective dielectric constants for s and p po- 
larizations 



'-eff 
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where c is the speed of light. 

In order to calculate the photonic band structure of the 
array of carbon nanotubes shown in figure la, we need 
to find solutions of the Maxwell equations with a tensor 
dielectric function. 



V X E = -/ioM dn/dt, 

V X H = +eo£(r,cj) d'E/dt, 



(4) 



that obey Bloch's theorem. On-shell methods like the 
one described in Ref. [8] are ideally suited to calculate 
the band structure of materials whose dielectric function 
depends on uj, as graphite does. In these methods, oj is 
fixed at first and hence e{uj) can be specified, allowing us 



to establish an eigenvalue equation for the calculation of 
k{uj). In particular this approach has a great advantage 
for studying EM waves in metals [|l5| where the dielectric 
function depends strongly on frequency. Very recently 
p6[ , an extension of this formalism has been developed 
which is able to work with tensor like dielectric functions 
as the ones present in eqs. (4). Details of this formalism 
can be found in Refs. [8] and [16]. 

With this methodology we calculate all the allowed 
Bloch waves of our system for a given photon energy uj. 
Then we can extract the dominant ones looking at the 
imaginary parts of their wavevectors. Finally, we can 
define the effective dielectric constants using Eq. (3). 

In its generalized version, MG approach is just the ap- 
plication of Clausius-Mossotti relations to describe the 
dielectric properties of composite materials in an effec- 
tive way. The basic assumption of this approach is that 
the polarizability of an individual object is modified by 
the presence of another ones only via their depolarizing 
contribution to the electric field acting on the object. For 
the case of arrays of cylinders in vacuum, a MG effective 
dielectric constant can be written as a function of the po- 
larizability of an individual tube a„{uj)^ and a depolar- 
ization factor Lct, which both depend on the polarization 
of the incident light, a: 



a.MG 



(^) = 1 + 



1 - fL^aaiio) ' 



(5) 



where / is the volume fraction occupied by the cylinders. 

For p-polarized light Lp — ^ and the polarizability of 
an individual plain cylinder is : 



where A 



ap(w) 



'eiiH + A' 



(6) 



Introducing this expression into 



Eq. (5) , the effective dielectric constant for p-polarization 
can be written as: 



pMG, , £||H + A + /(£||M-A) 
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For s-polarized light, Lg — and the polarizability of 
a plain tube is just the same as the polarizability of the 
homogeneous medium, as(w) = — 1, leading to a 

very simple formula for the effective dielectric constant: 



'-eff 



H=/£^M + (1-/). 



(8) 



If we wish to consider the general case of hollow cylin- 
ders, / must be replaced in eq. (8) by /' = /(I — 7^), 7 
being the ratio between the inner and outer radii of the 
tubes. 

Eq.(8) is the exact result for the effective dielectric con- 
stant for s-polarization. For this polarization there is no 
EM interaction between the tubes as MG approach as- 
sumes (Ls = 0) and the total polarizability is just the 



sum of the polarizabilities of the individual nanotubes. 
Because the electric field is directed along the tubes, the 
optical response of the system depends only on e±{Lu). 
Our numerical results for this polarization are in com- 
plete agreement with Eq.(8), giving additional support 
to our scheme. 

The simple relation (Eq.(8)) that holds for s- 
polarization can serve as an indirect way to estimate the 
volume fraction occupied by the nanotubes, /. The imag- 
inary part of J is simply /' multiplied by Ime . This is 
the reason why the experimental Ime^ f f Q has its max- 
imum at 4.6 eV just as Imej_ does Then, a rough 
estimation of / can be derived using the different heights 
of the peaks at 4.6 eV for the reported Ime^j^ and the 
experimental Ime_L . Taking the inner radius of the tubes 
to be 0.25 - 2 nm, we find / « 0.6 - 0.7. 

For p-polarized light, the EM interaction between the 
nanotubes is not negligible and our numerical method 
can give new insights in the analysis of this interaction. 
In figures 2 and 3 we show our numerical results for the 
effective dielectric constant, comparing them with the 
results of the MG-hke approach given by Eq.(7). 

With figure 2 our aim is to analyze the behaviour 
of these dielectric constants for low volume fractions, 
/ varying from / = 0.09 {d = 30nm) to / = 0.27 
{d = 17nm). As we can see from this figure, the good 
agreement between MG results and our method is re- 
markable for these values of /; only at low frequencies 
there are some small differences in Ime^^^ between both 
mean field theories. It is known that MG approach works 
very well in the dilute limit for dielectric media 0, and 
our numerical results show this is also the case for aligned 
carbon nanotubes. As we can conclude by looking at the 
two peaks present in line^jj:{uj), the optical response of 
the structure in the dilute limit is governed by the dipo- 
lar EM modes of an isolated carbon nanotube. Due to 
the anisotropy of the dielectric function of graphite, there 
are two dipolar modes, one located energetically at 5.3 
eV and the other at a higher energy, 6.3 eV. 

In figure 3 we study the case in which the EM interac- 
tion between the nanotubes becomes important, showing 
our results for the effective dielectric constant for dis- 
tances between the tubes smaller than three times their 
radii (/ > 0.3). In the same figure we show the results of 
Eq. (7) for this range of distances and the experimental 
data for p-polarization as reported in Ref. Q . We can see 
that the effect of the EM interaction is basically to shift 
to lower energies the low energy dipolar mode (from 5.3 
eV for isolated nanotubes to 4.8 eV for very close con- 
tact nanotubes), and conceal the dipolar mode at 6.3 eV. 
Again, MG approach seems to be a very good approxima- 
tion for e^j!j:{Lo) for intermediate distances {d > 11.5nm) 
revealing that, for this range of distances, the EM cou- 
pling between the tubes is basically a dipolar one as MG 
assumes. However, the inclusion of the full EM interac- 
tion between the tubes as our mean field theory does is 



necessary in order to describe properly the strong cou- 
pling limit (d < 11.5nm). The experimental imaginary 
part of the effective dielectric constant presents its max- 
imum just at 4.8 eV, and our results suggest that the 
position of this peak can only be explained if the nan- 
otubes are indeed in very close contact {d » 10.3nm) 
forming the films. As regards to the real part of s^ff, 
the shape of the theoretical curve for d = 10.3 nm closely 
matches the experimental data but a constant shift of 
approximately 1 is needed in order to obtain an excellent 
agreement between theory and experiment. 

In conclusion, we have presented an effective medium 
approach to analyze the optical properties of carbon nan- 
otube films. In contrast with usual effective medium 
methods, our numerical approach can treat the limit in 
which the EM interaction between the constituents of a 
complex system is important, as well as the dilute limit. 
We have also developed a Maxwell-Garnett approach to 
study this system. Using both mean field theories we 
have analyzed in detail the EM interaction between the 
nanotubes as a function of the volume fraction. MG ap- 
proach appears to be a good approximation for large and 
intermediate distances between the tubes. However, the 
inclusion of the full EM coupling between the nanotubes 
is essential in order to study closed-packed structures, ap- 
pearing in aligned carbon nanotube films ||^. The good 
agreement obtained between our numerical results and 
the experimental data demonstrates that our mean field 
theory can be a very useful tool to analyze spectroscopic 
properties of nanostructured materials. 
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Figure Captions 

Figure 1. a) Our model of a carbon nanotube filni: 
an array of cylinders, infinitely long in the z-direction, is 
arranged on a square lattice in the xy-plane. b) Unit cell 
of the structure, R = 5nm is the radius of the cylinders, 
d is the distance between the centers of the tubes and 
Ly is the thickness in the direction of propagation of the 
light that will be chosen equal to d. Also shown are the 
cylindrical unit vectors, 6 and f in order to describe the 
change of coordinates from cylindrical to Cartesians and 
the subsequent change in the tensor like dielectric func- 
tion associated to a generic point inside the tube (x,y) 
(see text). 

Figure 2. The real and imaginary parts of the effec- 
tive dielectric constant for p-polarized light, with d vary- 
ing from d = 30nm to d = 17nm. Full line: our numerical 
results. Dashed line: results using a MG approach (see 
eq.7). 

Figure 3. The real and imaginary parts of the effective 
dielectric constant for p-polarized light for different val- 
ues of d. (a) d = 10.3nm (b) d = lO.Tnm (c) d — ll.Snm 
(d) d = 12.5nm (e) d = 14nm. Full line: our numeri- 
cal results. Dashed line: results using a MG approach 
(see eq.7). The dots represent the experimental data as 
reported in 1^. For the sake of clarity, the curves for 
the real part are shifted by a quantity shown in the fig- 
ure. The experimental values of Ree^^j: are shifted by a 
constant value of -|-5 with respect to the reported ones. 
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